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In this work we generalize the Bochner criterion addressing the characteristic function, i.e., the 
Fourier transform, of the Glauber-Sudarshan phase-space function. For this purpose we extend 
the Bochner theorem by including derivatives of the characteristic function. The resulting neces¬ 
sary and sufficient nonclassicality criteria unify previously known moment-based criteria with those 
based on the characteristic function. For applications of the generalized nonclassicality probes, we 
provide direct sampling formulas for balanced homodyne detection. A squeezed vacuum state is 
experimentally realized and characterized with our method. This complete framework—theoretical 
unification, sampling approach, and experimental implementation—presents an efficient toolbox to 
characterize quantum states of light for applications in quantum technology. 

PACS numbers: 42.50.Ct, 42.50.Dv, 03.65.Wj, 02.30.Nw 


Introduction. Quantum physics differs in various 
ways from a classical description of nature. Exploit¬ 
ing these differences became a major research inter¬ 
est for the aim of developing quantum technologies. 
One principle scope of this held is the formulation of 
measurable conditions, which are fulhlled for classical 
systems but may be violated for nonclassical systems. 

In quantum optics, a fundamental approach for sep¬ 
arating classical light helds from quantum ones is 
dehned in terms of the Glauber-Sudarshan P func¬ 
tion mm- If a harmonic oscillator system can be 
solely described by coherent states and classical statis¬ 
tics, then its P function is a classical probability den¬ 
sity. Whenever the P function does not exhibit such 
classical behavior, it identihes a nonclassical quantum 
state [5]. In general, the P function can represent a 
quasiprobability. For example, squeezed light cannot 
be characterized by a classical P function, as has been 
experimentally conhrmed 015]. 

In order to certify quantum features of states, a 
number of nonclassicality criteria have been formu¬ 
lated. Here, let us focus on two hierarchies of criteria 
which have been successfully applied in theory and ex¬ 
periment to determine the quantum character of light. 
Both of them prove the nonpositivity of the P func¬ 
tion in terms of experimentally accessible quantities. 

The first hierarchy is based on Bochner’s theorem 
and the characteristic function (CF), being the Fourier 
transform of the P function [ail]. Bochner’s theo¬ 
rem expresses necessary and sufficient conditions for 
a function to be the Fourier transform of a classical 
probability density [5]. It presupposes only the fact 
of the function to be defined for all arguments, so 
that it is nonlocal in phase space. The CF can be 
directly sampled from balanced homodyne detection 
(BHD) data, as shown in experiments [5fH2). uncover¬ 
ing nonclassical phenomena through violations of the 
Bochner conditions. From the hierarchy of the latter, 
only the second-order condition has a clear geometri¬ 
cal interpretation. Higher-order conditions, however, 
fail to exhibit such a geometric meaning, since they 


correlate a higher number of points in phase space. 

The second hierarchy is formulated in terms of the 
so-called matrix of moments (MOM) |T3l |14|, con¬ 
taining the statistical moments of the P function. 
Some second-order minors of MOM are known to 
identify fundamental quantum effects, such as sub- 
Poisson [inj or squeezed light [TB]. The MOM, how¬ 
ever, includes more general quantum effects in terms 
of higher-order moments (cf., e.g., [TTHSO]). Generaf- 
izations of this method even allow one to identify en¬ 
tanglement inum and space-time-dependent quan¬ 
tum correlations [53|. 

In the present contribution we combine the advan¬ 
tages of the GF and the MOM of the P function re¬ 
sulting in a generalization of Bochner’s theorem. This 
leads to a hierarchy of nonclassicality conditions em¬ 
bedding CF and MOM nonclassicality probes. Be¬ 
yond the established cases, we obtain previously un¬ 
known conditions, which are shown to outperform the 
previously known criteria for specific mixed quantum 
states. In order to apply the proposed technique, 
we derive sampling formulas to reconstruct deriva¬ 
tives of the CF from BHD data. The path—starting 
with the theoretical treatment of unified nonclassical¬ 
ity probes—is finalized with our experimental imple¬ 
mentation of the criteria for a squeezed state. Even 
though such states have a demanding sampling be¬ 
havior, due to the exponential growth of their CF, 
the nonclassicality of the measured state is confirmed 
with a high significance. 

Identification of nonclassical CPs. Any single¬ 
mode quantum state of light can be given in Glauber- 
Sudarshan representation 00 , 

p = j d^aP{a)\a){a\, (1) 

where P is a quasiprobability distribution. On the 
one hand, if P is a classical probability density, the 
state p has a classical counterpart. On the other hand, 
nonclassical states are characterized through negativ¬ 
ities in the quasiprobability distribution, such as Fock 
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or squeezed states. The P function is often strongly 
singular [5] and, hence, not always experimentally ac¬ 
cessible. 

Even if the P function is a highly nontrivial distri¬ 
bution, its CF, 

$(/?) = = J d^a (2) 

is a well-behaved function. Here : • : denotes the nor¬ 
mal ordering prescription, and d and are the annihi¬ 
lation and creation operators, respectively. Bochner’s 
conditions [8]— for $(/3) to be the CF of a classical 
probability distribution—were reformulated for quan¬ 
tum optics in [J. Namely, $(/3) is the CF of a classical 
state if the following conditions are fulfilled: (i) nor¬ 
malization, $(0) = 1; (ii) hermiticity, <&(—/?) = $*(/?) 
for all /3 S C; and (iii) positive semidefiniteness: for 
any positive integer N and arbitrary complex numbers 
Pi,...,I5n holds 

<i> = [d>(/3i-/3,)]5=i>0, (3) 

i.e., the matrix $ is always positive semidefinite. Since 
conditions (i) and (ii) are valid for any state, the vio¬ 
lation of condition (iii) discerns classical and nonclas- 
sical states. 

Employing Sylvester’s theorem, the non-negativity 
in (§ can be written in terms of minors. For example, 
the second-order nonclassicality condition {N = 2) 
reads as 

det$ = 1 - |$(/3)P < 0, (4) 

having a clear geometric interpretation. Whenever the 
absolute value of the CF exceeds 1, the state is non- 
classical [B]. Beyond this second-order minor, there is 
no clear geometric view of the higher-order conditions. 

Generalizing Bochner’s theorem. For the desired 
generalization, let us recall that for any classical state 
and any function f{a) holds 

y d^aP(a)|/(a)p > 0. (5) 

The relation of the inner product of two functions g{a) 
and h{a) to the inner product of their two-dimen¬ 
sional Fourier transformed functions g{f3) and h{/3) is 
given from Parseval’s theorem: J d?ag*{a)h{a) = 

J d‘^f3g*{l3)h{l3). Hence, the non-negativity condi¬ 
tion can be written as 

yd2/3$(/3)Aut[/](/3)>0, (6) 


Now, the classicality condition reads as 


N 

E 

*j=i 




rii -\-mi 


$(/?) 


/3=/3.-ft 


> 0 , 


( 8 ) 

(see the Supplemental Material [24] for details). 

The definition of a generalized Bochner ma¬ 
trix (GBM), 


5$ = 


(-ir 


ni+miQrii+mj j^rij+m 


dn 
fl 


‘^ci>(/3) 




1 N 




(9) 

and the vector / = (/i,..., /at)"'’ leads to the gener¬ 
alized version of Bochner’s theorem. Before formulat¬ 
ing the details, let us comment on some properties of 
5$. First, the GBM is Hermitian [23]. Second, in¬ 
equality Q may be formulated in the compact form 
> 0. Third, the elements of GBM can be also 


given as 








( 10 ) 


Theorem. For any classical state it holds that the 
GBM is positive semidefinite, i.e., 9$ >0. ■ 

Properties of the generalized theorem. As men¬ 
tioned earlier, a minor representation for the gener¬ 
alized Bochner theorem can be formulated. Namely, 
a state is nonclassical, if there exists a positive in¬ 
teger N, non-negative integers n = {ni,... 
and m = (mi,..., mAr)"^, and complex numbers /3 = 
(/?!,..., /SAf)"*" such that 

det 5$ < 0. (11) 


As an example, let us construct a nonobvious non¬ 
classicality criterion. The parameters N = 3, rh = 
(0,1,0)"^, n = (0,0,1)"'’, and P = (/3,0, for /3 G C 
produce the GBM: 


/ $(0) dp^ld) 5;3.$(0) \ 

-dp.^-13) -dpdp.^O) - 9 |.$(-/ 3 ) . 

V -9^$(0) -a2$(/3) -dgdp.^0)J 

( 12 ) 


This previously unknown nonclassicality criterion 0 
for the GBM (12) combines moments and the GF of a 
nonclassical state. That is, for /3 = 0, the minor corre¬ 
sponds to the quadrature squeezing condition mm 


using the Fourier transform of the function |/(a)p 
which is given by the autocorrelation function 
Aut[/](/?) = f d^jf*(j)f('y + /3). Choosing arbitrary 
non-negative integers irii and as well as complex 
numbers Pi and /i, we can define a function 

N 

m = ( 7 ) 


det{d^) 




,)]2:)(:[Ai(^^a,,)]2:)<0, 

(13) 

using the relation (10). The general relation between 
CF and MOM nonclassicality conditions will be de¬ 
rived in the following by applying the above theorem. 

First, we may compare the GBM, Eq. (|^, with 
the corresponding MOM and Bochner matrices of the 
same dimensionality N. A Bochner matrix $, defined 
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in Eq. (§. corresponds to a GBM 9$ in Eq. if we 
set the integers as mi = ■ ■ ■ = itin = ni = ■ ■ ■ = un = 
0. Thus, the matrix elements of the GBM reduce to 




/3=/3.-A 


= $(A-/3,). (14) 


Since constraints (i) and (ii) of a GF are already ful¬ 
filled, this special case of the non-negativity of the 
GBM is equivalent to condition (iii) in inequality (§. 
Gonsequently, 9$ > 0 implies $ > 0, which is al¬ 
ready necessary and sufficient to characterize a classi¬ 
cal state [7j. Note that the Bochner criterion is non¬ 
local in phase space, since it correlates distant points 
Pi and Pj (* 7^ j). 

Second, a MOM corresponds to a GBM, if we choose 
/3i = • • • = /Sat = 0 in (10), since the matrix elements 
simplify to 5$ = [(-1)"^+". 

Using the representation of the non-negativity in 
Eq. (Isl), we obtain 


N 

*.i=i 


(15) 


with f'j, = (—l)”'=/fc for k = 1,N. In this special 
case, the generalized Bochner theorem coincides with 
a MOM criterion for annihilation and creation oper¬ 
ators [14] . Since the MOM is represented through 
derivatives of the CF at the origin, /3 = 0, they are 
local characteristics in phase space. 

The previously independent approaches of the GF 
and the MOM have been unihed and generalized by 
our method. Nonlocal Bochner and local MOM non- 
classicality probes occur simultaneously in the non- 
classicality theorem based on the GBM. In the fol¬ 
lowing, we provide an example of a particular GBM 
inferring quantumness of a state which is inaccessible 
via the corresponding MOM and Bochner approaches. 

Let us consider the situation, when the condition (|^ 
does not verify the nonclassicality of a given state. 
Additionally, the MOM minor 


det ( , 1 > 0, 

(aT) (:aTa:), 


(16) 


is non-negative for any state, as it represents the mean 
number of incoherent photons M- The second-order 
minor resulting from the GBM and the choice mi = 
m 2 = ni = 0, 712 = 1 and P = Pi — P 2 reads as: 


det 9$ = det 


1 —(:ae 


b) 


:a’^a:) 


= (:afa:)- = ([Aif][Ai]), 


(17) 


with A = “. Now, we construct a mixed 

state, Pniix, being a mixture of a thermal state with 
three- and four-photon-added thermal states [2HH27] . 
Its CF is 


<i>n.ix(/?) = + P3L3[(1 + h3)|/3p]e-«^l^l' 

+ P4T4[(l + h4)|/3p]e-^-l^l', (18) 



lm(P) -10 -10 Re(P) 


FIG. 1: (Color online) Left half: The scaled (3.06x) 
minor IFt for the state ( |18[ ) is shown for the parame¬ 
ters po = 0.944, ho = 0.1, P 3 = 0.03, ha = 0.12, p 4 = 
0.026, h 4 = 0.182. The displayed negativity at |/3| « 5.8 
of this determinant verifies the nonclassicality of the state 
under study. Right half: The lowest-order criterion Q 
does not verify the nonclassical nature of the state. 


with probabilities pi, mean numbers of thermal pho¬ 
tons fii, and the Laguerre polynomials Lpz) (i = 
0,3,4). The visualization of the minor 0 for this 
state is given in Fig.[2 It can be seen that negativities 
infer the quantumness of this state, which can neither 
be observed through the nonclassicality condition ([^ 
nor the MOM minor (16) (see [24] for details). 


Sampling the CF and its derivatives. Properties 
of quantum states of light may be directly sampled 
from BHD [28l [29] . In this setup one measures the 
statistics of the field quadrature operator, x{(p) = 
e*‘^a -I- e“*‘^a1', resulting in M data points of quadra¬ 
ture values [xj{p)]^i for a fixed phase p. The CF (H 
is readily obtained from the set of data. 


$(/3) « (19) 

1=1 


where arg/3 = 7r/2 — (cf. [5HTT]i. For sampling 
the derivatives in an analogous way, the knowledge of 
two noncommuting operators, [d, fll] = 1, is required. 
Hence, another ansatz for reconstructing the GBM 
elements has to be developed. 

In order to apply our generalized Bochner condi¬ 
tions to experimental data, we need to reconstruct 
the characteristic function and its derivatives. An al¬ 
ternative way to reconstruct the CF is based on its 
Taylor expansion. 


4>(/3) = (: :) = £ ( 20 ) 

fc ,/—0 

The normally ordered moments can be represented 
with quadrature distributions p(x, ip) and Hermite 
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polynomials Hn{z) as 


(a^‘a‘ 


TT 


e^{k-i)vf.ii\ 


V¥^{k + l)\ \V2 


( 21 ) 


Hk+i I ^ 


The double series (|20|) can be evaluated as 

p{x,(p) 


pTZ p-\-00 

= dip dx 
Jo J — oo 


7r(7 + 7*) 

X + 7’e 


( 22 ) 

—X'Y* —7*^/2 


with 7 = (see 


for more technical details). 


By differentiation of (22 1 the following pattern func¬ 
tion emerges: 


^ r!(-l)«+'=i+'=3(9 + /ci)!2-'=3/2 

,(a:,7)- 2^ fci!fc2!A:3!(7-h7*)9+'=i+i ‘ 

fcl -|-fc2 
-\-k3—r 

X exp(a:7 - /2)Hk^{j/V2), 


,k2 


(23) 


and = 0 for q<0 or r’<0. Now, the derivatives of 
the CF can be directly sampled as 


(24) 

i=i 

X [mD ^-1 {x ,, /3e*^0 + 


with the need of a uniform distribution of the mea¬ 
sured data {xj,ipj)Y^i with respect to the phase ip. 

Experimental implementation. We generated the 
squeezed field through a hemilithic, standing wave, 
nonlinear cavity (see Fig.j^ serving as an optical para¬ 
metric amplifier (OPA). As nonlinear medium 
we used an 8 -mm-long, 7% magnesium-oxide-doped 
lithium niobate ( 7 %Mg 0 :LiNb 03 ) crystal. The OPA 
is pumped with a strong laser beam (290 mW) at 
532 nm, resulting in a classical gain of 3.3, which 
yields —4.13dB squeezing and 6.11 dB anti-squeezing 
at 1064 nm. A BHD was implemented [HU |33] with 
98% visibility and a quantum efficiency of 90%. Along 
with additional losses in the squeezed field due to 
losses in optical components and an escape efficiency 
smaller than unity, this results in a overall efficiency 
of 77%. The optical phase of the signal, relative to 
the local oscillator (TO), was controlled by a mirror 
mounted on a piezoactuator. A continuous variation 
of the optical phase ensures an equally distributed 
phase of the quadrature data points {xj,ipj)^i [5]. 


Using the sampling formula (24), we reconstructed 
the determinant of the GBM (12). It is important to 
mention that the CF of the squeezed state is diverging 
in one direction of /?, which causes high sampling er¬ 
rors. This fact makes such a state an optimal test for 
demonstrating the applicability of our sampling ap¬ 
proach under difficult premises. Figure shows the 






FIG. 2: (Color online) Experimental scheme for the gen¬ 
eration and measurement of a squeezed state. The state 
was generated in an OPA with an 7 %Mg 0 :LiNb 03 crystal, 
which was pumped with a 290 mW laser beam at 532 nm 
generating —4.13dB squeezing at 1064 nm. A BHD (LO 
power 1.23 mW) is set up for the detection. The phase of 
the LO was altered continuously by applying a triangular 
alternating voltage onto the piezo that is responsible for 
the mirror position. 


determinant det(9$) of the GBM, Eq. (12), together 
with its signed significance. 


E[det(54>)] = 


det(94>) 

(T[det(5$)]’ 


(25) 


where cr[det((9$)] denotes the standard deviation. 
Without an additional search for optimal phases, the 


squeezing condition (13) is fulfilled as det(9<i))|^=o = 
—0.469 ± 0.007. At this point the highest significance 
of a negative value of above 70 standard deviations is 
reached. From the decay of the significances, it can 
be observed that the squeezed and anti-squeezed axes 
are—up to a slight rotation—the imaginary and real 
axes, respectively. 

Conclusions. We derived necessary and sufficient 
nonclassicality probes by generalizing Bochner’s the¬ 
orem. This generalization has been proven to em¬ 
bed two hierarchies of previously known nonclassical¬ 
ity criteria. Namely, the conditions resulting from 
the original Bochner’s theorem and the matrix of 
moments can be reproduced by suitable parameter 
choice. Our generalized Bochner theorem makes use 
of the advantages of the original Bochner theorem, 
being nonlocal characteristics in phase space, and of 
the matrix of moments, directly yielding quantum fea¬ 
tures such as squeezing. 

We explicitly constructed a nonclassical state, for 
which specific second-order moments and Bochner cri¬ 
teria fail, but the corresponding generalized second- 
order criterion truly visualizes its quantumness. Pat¬ 
tern functions have been derived to sample our cri¬ 
teria from balanced homodyne detection data. For a 
direct demonstration, a squeezed state—having a de¬ 
manding sampling behavior—has been prepared and 
measured. We could successfully demonstrate that 
our approach certifies the nonclassicality of this state 
with high significance. 

Thus, we explored the entire process from formulat¬ 
ing a unified technique to uncover quantumness un¬ 
til its successful experimental implementation within 
this work. Moreover, it is straightforward to gener¬ 
alize our method to multimode scenarios by applying 
the present techniques. The here established method 
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FIG. 3: (Color online) The top plot depicts the experimen¬ 
tal implementation of the nonclassicality condition for the 
determinant of (121, det((?4>) < 0. Negative values demon¬ 
strate that the GBM criterion is capable to prove the non¬ 
classicality of the generated state. The bottom plot shows 
the signed significance E[det(9<I>)] = det(9<E')/a[det(9<I>)] 
up to —70 standard deviations at the origin and -|-60 in 
side peaks. The region with |E| < 5 is shown white. The 
error CT[det(94?)] is obtained via a linear error propagation 
of the GBM elements. 


allows one to infer quantum effects of previously un¬ 
known structures in phase space beyond local prop¬ 
erties, such as higher-order moments, and nonlocal 
properties of the characteristic function, as addressed 
by the original form of the Bochner theorem. Even¬ 
tually, this unified characterization of quantum fea¬ 
tures will allow one to assist quantum technologies 
with properly designed quantum states of light. 
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Section I includes details on the generalized Bochner matrix. In Section II we discuss the physical 
meaning of some low-order minors. We also highlight properties of special kinds of states. In 
Section III the detailed derivation of sampling formulas and the proof of pattern functions to be 
analytical is given. Additional information about the experimentally realized squeezed state and 
the condition used to prove its quantumness is shown in Section IV. 

We use the abbreviations of the main part: generalized Bochner matrix (GBM), characteristic 
function (CF), and matrix of moments (MoM). 


I. THE GENERALIZED BOCHNER MATRIX 


In this section the GBM is derived and its properties are studied. The special cases of the Bochner matrix and the 
MoM are retrieved. 

The autocorrelation function of a function /(/?) = ~ Pi) reads as 


Aut[/](/3) 



( 1 ) 


Using integration by parts we get 


JV 

Aut[/](/) = ^ J + 7-/3/ 


ui=i 

N 


= E +Pi - Pi)- 

ij=l 


( 2 ) 


The evaluation of the inner product / of this autocorrelation function with the CF $(/3) yields 

1= J d^p^p)Ant[f]{P)= J d^PHP)d;^^"^^dp+^^6{p + pj-Pi) 

i,j=i- 

N N 


(3) 


i,f = l 


bi=i 


P=Pi-Pj 


Defining the vector / = (/i,..., /jv)"^ and the GBM 

'9$ii ... 


5$ = 

I 

••• d^N,N) 

we can write the inner product I as 


= 1 


3=ft- + 


bi=i 




(4) 


( 5 ) 
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The elements of the GBM can be also computed as 




ni+rij /.;:,fni+m,j ~mi+nj l(}i-0j]a^-HSi-Pj]" a.\ 


/3=/3i-ft 




The GBM is Hermitian, because for all i,j holds 

d^* . =(— 

For the special case n* = mi = 0 for all i, we get the corresponding Bochner matrix 


0$ = 


/ $( 0 ) 

Hh-Pi) 


HPi- 

m 



Note that trp = 1 = f (PaP{a) = ^*(0) for any quantum state p. 

For relating the GBM to the corresponding MoM, one may generate a list of pairs of natural numbers; 


0 \ l\ / 0 

oj’loj’lu’Vo 


0 


0 

n- 1 y ’ V 0 


n — 1 
1 


1 

n — 1 


0 \ / n + 1 


The upper and lower vector components yields an ordering as 

(ni, n2,...) = (0,1,0, 2 , 1,0,...) and (mi, m 2 ,...) = (0,0,1,0,1, 2 , ...). 
For /3i = 0 for all i, the GBM corresponds to the MoM [1], 


d^ = D 


/ 1 (a) (at) •••\ 

(at) (ata) (at^) 


(o) (a^) (ota' 

V ; : -.y 

up to the unitary transformation D = diag ((~l)”0ieN\{0} ~ ~ D~^. 


D, 


II. LOW-ORDER MINORS AND SOME EXAMPLES 


A. Interpretation of low-order minors 


Second and third order minors can be considered as the variance or cross-correlation, respectively. This 

{Ax = x-{x)y. 


det ={A*A) - {A*){A) = (AA*AA), 

( 1 (^) {B) \ 

det (A*) {A*A) {A*B) ={AA*AA){AB*AB) - {AA*AB){AB*AA). 
\{B*) {B*A) {B*B)j 

Second order nonclassicality probes of the Bochner criterion or a specific MoM criterion are 

=1-|<1>(/?)P<0, 

^®^((at) =([Aa]t[Aa])<0. 


( 6 ) 

( 7 ) 

( 8 ) 

( 9 ) 

( 10 ) 

means 

(11) 

( 12 ) 

(13) 


(14) 
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B. Multi-photon-added-thermal-states. 


Note that states with a Fock diagonal statistical operator p have a rotationally symmetric CF, $(/?) = <F(|/3|). One 
class of such kind of states are multi-photon-added-thermal-states [2]. The P function of a m photon-added thermal 
state (mPATS) is 


P{a) 


(-ir 




1 

-^e 

Tin 


5 


(15) 


where n is the mean photon number of the thermal background and is the mth Laguerre polynomial. The latter 
can be expressed by the Rodrigues formula: Lm{z) = The Fourier transform of the P function of 

the mPATS PmPATs{n) is 

cl>(/3) = L„[(l + n)|/3|2].e-^l'3|^ (16) 


The single-photon-added-thermal-state (SPATS or, in our notion, IPATS) has been experimentally realized in [3] 
and analyzed in [4]. We use a statistical mixture of a thermal state (OPATS), a 3PATS, and a 4PATS to show the 
advantage of the GBM approach. That is 


Pmix = PoPoPATs(?^o) +P3P3PATs(»^3) + P4P4PATs(»^4), 


(17) 


with the characteristic function: 

^mix(^)=Poe-^«l^l'+P3L3[(l+n3)|/?P]e-«^l'5|'+P4L4[(l + n4)|/3|2]e-’^-l^l'. (18) 


The chosen parameters are the probabilities po = 0.944, ps = 0.03,p 4 = 0.026 and mean numbers of thermal photons 
ho = 0.1, he = 0.12, hg = 0.182. The large contribution of the classical thermal state, po ^ P 3 ,P 4 , does not allow 
to proof the nonclassicality of this state with Bochner condition in Eq. (13); see dashed line in Fig. 1. The MoM 
criterion in (14) reads as ((a — (a))^(d— (a))), which cannot be negative since X^X > 0 for all X. However, by using 
the corresponding determinant of the GBM in the main body. 


det 9$ = det 


{:dpe- 


1 , 


•ae^ 


we can prove the state (17) to be nonclassical; see the solid line in Fig. 1. 

det(^O) 



1/31 


(19) 


FIG. 1: Solid line shows the radial course of the GBM determinant (19). It proves nonclassicality for |/3| 5.8. Both, the 

Bochner determinant in (13) (dashed line) and the MoM determinant in (14) (the value of solid line at \j3\ = 0), are non-negative 
for this state. 


III. SAMPLING FORMULAS 

We will rewrite the GF as an average of the quadrature distribution for any phase. This renders it possible to 
formulate a sampling approach for the GBM. 
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Under the assumption that all moments exist, we can expand the CF as 

m = :) = £ 

fc,Z=0 

The normally ordered moments {a)^a}) can be reconstructed from the quadrature distribution p[x,<f) [5]; 

, r+oo f x \ 


pTT P-rOC 

= / d(p dxp{x,(fi)- 
Jo J-oo z 


' ' Jo ^J-oo ’"VV^(m + n)! “■ 

where is the fcth Hermite polynomial. Inserting Eq. (21) into Eq. (20) yields: 


Hm+ra 


x \ 


^TT /*+00 

^(/^) — 

Jo J-oo 




Jo ^ J-oo - V 2 ^(fe + ^)! --‘Vvsy 

Now, the latter sum may be rewritten (7 = /le**^): 

“ ^ X - (_^)P ( x\^\ 7 r (- 7 )" „ / ^ ^ 1 - (- 7 / 7 *)^+' 


7 t-7 1 c - ' u ^ 1 1 7 tj / j(_ nr -j- -’ntizizIh / j(_ irAz2zj_i__ to'?'! 

^^(fc + 0! l-(-7/r) ^ 

The remaining series can be directly calculated with the generating function of Hermite polynomials [ 6 ]: 

^ fU 

^-H„(a:) =e2"*-A (24) 

n=0 

Inserting this relation into the CF, it can be expressed with quadrature distribution p(x, (f>) and 7 = I3e^^ as 

m= dp dxP^J^ --. (25) 


.7 0 J — 00 


7 + 7 * 


For the proof of the integrand to be analytical, we consider the expansion (7 = 7 ij + iji)- 

.yga:7-7^/2 _|_ ,y*g-a:7*-7*^/2 ga;7-7^/2 _ g-a:7*-7*^/2 

/(O') = -UTT 7 - = A-—-- 


7 + 7* 


“ — e 
7 + 7 * 


—d g»^7J-7g/2+7?/2 ^^’^^[0'-R(^ W)] _ Q ia:77-7^/2+7f/2/ _ ■ 7 , hnix ill)] 

IR ■ " ^ ( 2 * + l)' ■ 

Note that the remaining series converges absolutely: 

00 I . . N| 2 fc 00 I . , .| 2 fc 

^ ( 2 fc + l)! - ^- m- -= ■ 


cxD I / . x| 2 fc 00 I / . x| 2 fc 

7i^(2:-^7/) ^^\lR{x-i'yi)\ 

L ^ L-rn-= i7H(2: - *7r)|. 


The power series of the integrand converges for all 7 , so the integrand is an analytical function. 

Let us now derive such a representation for derivatives of the CF. Using Eq. (25) with 7 = we need to replace: 




Now, arbitrary derivatives can be computed: 


d^d]^M)= [ dp r , 

^ ' Jo J-oo ^ ^ ^ [7 + 7* 


77-7^72 _|_ O' g-a: 7 ’‘- 7*'®/2 


f\p 

j'+oo 

da; ghm-n)v 

Jo J 

— 00 

TT 

f dp 

p+oo 

da; fp(m-n)^p 

Jo J 

f 

— 00 

IX 


^_1 (-l)”n! ^,r 7 - 7 V 2 


(^ _j_ ^*^ra+l 


7 + 7* 


' 7 Z^_|_^*Xn+l 


„_1 7 ._ 7-/2 , (-l)"^m! 




5 


Using the definition of Hermite polynomials, H„(2/\/2) = (—^ i 


we define the pattern functions: 




E 

ki+k2+k3.=r 


H (_i).+fc.+fc3(g + fc^)!2-fe3/2 

ki\k2\k3\ (^ + ^*)9+fei+i 



(31) 


which has to be set to zero -DJ = 0 for g < 0 or r < 0. After evaluating all terms of derivatives the derivatives of the 
CF reads as 

dfd^0,^{P)= [ dip r + 

Jo J—oo ^ 

+ nDl-\-x,l3*e-^n + P*e-^‘^D^i-x,l3*e-^n] ■ (32) 

Using the standard sampling approach for quadrature-phase data points {xj,pj)jLi, we get the estimate 

1 “ 
i=i 


IV. EXPERIMENTAL RESULTS 



FIG. 2: Using Eq. (33) for m = n = 0, the sampled CF, Re[<I>(/3)], is shown. The imaginary part is zero within the sampling 
error, Im[<l>(/3)] 0. It is important to stress that the behavior along the imaginary axis (Re(/3) 0) has an inverse Gaussian 

slope. 


The sampled CF of the measured state is shown in Fig. 2. This squeezed vacuum state is theoretically described by 

$,,(^) = e-R'=(/3-^CV„/2+Im(/3e-)V./2+|/3|V2 ^ jg, (34) 

where 14 and 14 are squeezed and anti-squeezed variances, respectively, and p describes a rotation of the principal- 
axes (squeezed and anti-squeezed quadratures) in phase-space. The experimentally obtained variances are 14 = 4.086 
(or 6.11 dB anti-squeezing) and 14 = 0.386 (or —4.13dB squeezing). Note that for vacuum holds 14 = 14 = 1. For 
N = 3 and the parameters (mi,m 2 , m 3 ) = (0,1,0), (ni,n 2 ,n 3 ) = (0,0,1), and {Pi,p 2 ,P'i) = (/?,0,/3), the GBM in (4) 
reads as 

/ $(0) dp^{p) 5^.$(0) \ 

-dp.^{-p) -dpdp.m -dlM-P) 

\ -dp^iQ) -dpdp^m) 


= 


( 35 ) 







6 


Using (10) the determinant for /3 = 0 can be rewritten as; 


det(9$) 


/1 

(at) 

(a) \ 

= det (a) 

(at a) 

(a') 

/3=0 \ / k|\ 

(dt2) 

(ata)/ 


-(:[Ai((pn 




(36) 


being known from quadrature squeezing condition [1]. The last equality holds only for Gaussian states, which are 
fully characterized by two variances Va and I 4 . The sampled minor at the origin is det(cl$)|,g=o = —0.469 ± 0.007. 
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